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Abstract 

The long-time evolution of the KdV-type solitons propagating in ferro- 
magnetic materials is considered trough a multi-time formalism, it is governed 
by all equations of the KdV Hierarchy. The scaling coefficients of the higher 
order time variables are explicitly computed in terms of the physical parame- 
ters, showing that the KdV asymptotic is valid only when the angle between 
the propagation direction and the external magnetic field is large enough. 
The one-soliton solution of the KdV hierarchy is written down in terms of 
the physical parameters. A maximum value of the soliton parameter is deter- 
mined, above which the perturbative approach is not valid. Below this value, 
the KdV soliton conserves its properties during an infinite propagation time. 
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1 Introduction 



1.1 KdV-type solitons in ferromagnets 

Electromagnetic wave propagation in ferromagnetic media is intrinsically 
nonlinear. It is therefore the matter of intensive research in the theoret- 
ical physics of the nonlinear waves. In the frame of the Maxwell-Landau 
model, analytical expressions describing solitary wave propagation out from 
any slowly envelope or long- wave approximation have been found pQ. These 
waves have also been studied numerically Envelope solitons have been 
studied from several theoretical approaches jH [3] • There are many experi- 
ments regarding magnetostatic waves in thin films 00 EI- Long- wave type 
approximations allow to describe some features related to relativistic domain 
wall propagation [HllHj, but have also brought forward the existence of an- 
other type of wave, described by the Korteweg-de Vries (KdV) equation [TU] . 
It has been shown that such a wave can be emitted by a transverse instability 
of the relativistic domain wall [3 E|. The interaction between the two types 
of wave has also been studied [TU] . 

The KdV model is obviously a rough approximation. In ^UJ where it has 
been first derived in this frame, anisotropy, damping and inhomogeneous ex- 
change were neglected. Second it assumes that the wave depends on a single 
spatial coordinate (plane wave), and that the amplitude is weak enough, the 
wave length and the propagation distance large enough, so that the first or- 
der of the KdV approximation can be retained. A study taking into account 
three space dimensions, damping, and inhomogeneous exchange is published 
independently [11] . But the weakly nonlinear approximation itself may neces- 
sitate higher order corrections. The latter are independent from the former 
ones. Indeed, the wave is intrinsically nonlinear, and the weakly nonlinear 
approximation is forced by the introduction of a static field, to which the 
wave field can be compared. Even in the roughest approximation, the ra- 
tio between the two fields can become rather close to one. A derivation of 
the equations describing the evolution of the higher order terms has been 
derived using a multi-time formalism [T2] . It allows to prove that a formal 
asymptotic expansion exists up to any order, with all its terms bounded [J3] , 
which is a first step in the mathematical justification of the convergence of 
the expansion. Following the idea by Kraenkel et al. |14[ [15] . the multi- 
times expansion for KdV uses the KdV Hierarchy. The evolution of the main 
term in the expansion relative to each higher order time variable is given by 
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the corresponding equation in the KdV Hierarchy. Regarding the main term 
only, all information about the particular physical situation considered is 
contained in scaling coefficients of the time variables. These coefficients can 
be computed. The aim of this paper is to give the value of these quantities, 
and to draw physical consequences from them. It is organized as follows: in 
section 2 we describe the perturbative scheme in the multi-time formalism. 
In section 3 we compute explicitly the time scaling coefficients. Conclusions 
can be drawn on the validity of the pertubative scheme considered as an 
asymptotic expansion, i.e. for a fixed number of terms, when the perturba- 
tive parameter becomes small enough. In section 4, we give the expression 
of the one-soliton solution of the complete KdV hierarchy. This gives infor- 
mation about the validity of the perturbative scheme considered as a series 
expansion, i.e. for a fixed value of the perturbative parameter and an infinite 
number of terms. 



2 The multi-time formalism 
2.1 The KdV mode 

The evolution of the magnetization density M in a magnetic field H is de- 
scribed by the Landau equation 

d t M = -7// M AH eff , (1) 

where 7 is the gyromagnetic ratio (7 > 0), and /xo the magnetic permeability 

— * 

in vacuum. The effective field i? e // contains several terms giving account 
for the inhomogeneous exchange interaction, the effects of finite size and the 
anisotropy. Here we use the basic approximation: H e ff = H. Damping is 
also neglected. 

— * 

The evolution of the magnetic field H is described by the Maxwell equa- 
tions. We assume that, regarding its dielectric properties, the material is 
perfectly linear and isotropic, and we denote by c the light velocity based 
on its dielectric constant e, i.e. c= 1/y/e/io. The Maxwell equations reduce 
then to 

- V (V • Hj + AH = ^8 2 t (if + Af) . (2) 

We replace below H, M and t by the normalized quantities 7/i if/c, 7/UoM /c 
and ct. The constants 7/x and c take then the value 1. 
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The 'long wave' limit of a wave with negative helicity is considered. We 
introduce a small parameter e, such that 1/e measures the length of the 
solitary wave and e 2 its amplitude. The magnetic field is expanded as 

H = H + e 2 H 2 + ■■■, (3) 

and M in an analogous way. Using the slow variables 

n= eh, w 

it is shown first that this wave propagates at the velocity 



V= ^/(a + sin 2 #)/(a + l), (5) 

where 9 is the angle between the propagation direction and the applied field, 
and a = Hq/Mq the ratio from the latter to the saturation magnetization. 
Second it is shown that the propagation of this type of 'long waves' is gov- 
erned by the KdV equation [TO] 

d Tl (f2 + <?y?2<%y? 2 + rdfip 2 = 0, (6) 
where q and r are real constants given by 

_ 3 cos 2 9 sin 2 6y/l + a 

q ~ 9 '/ , • 2^3/2 ' ' ' 

A (a + sin 2 9) 

and 



— 1 cos 4 9\J a + sin 2 9 
r= 2^ sm 2 9(l + a) 7/2 ' U 

if 2 is the wave amplitude, related to the main component H2 and M2 of the 
wave magnetic field and the magnetization density through 

H 2 = ip 2 hi and M 2 = ^2^1, (9) 

where h\ and rhi are polarization vectors defined by 

sin 6 cos 8 
a+sin^ 6 

h 1 = m(l + a)sm9 | 1 | , (10) 




and 

m (1 + a) sin 6* cos # . . 

mi = . 2 cos# . (11) 

a + sin v 

(We use the normalization of piJJ ^] , introduced for computational conve- 
nience) . 

2.2 Higher order terms 

Going further in the resolution of the perturbative scheme, it is seen that the 
field component of order j (j > 2) writes 

Hj = (fjhi + h° (tp 2 , <pa,---, <pj-i) , (12) 

where hj is an explicit functional of the lower order amplitudes ip2, <f3, up 
to <fj-i, hi the polarization vector given above by (JTUJ), and tpj is an higher 
order amplitude, (fj satisfies a linearized KdV equation of the form 

d n <fj + qd ( (^j ) + rdf^j = 

Z,j{tp2,<P3,---<Pj-l), 

where the right-hand-side (rhs) member Hj is an explicit functional of the 
lower order amplitudes <£> 2 , ^3; U P to (pj-i- The parity and homogeneity 
properties of the expansion allow to prove that half of these equations admit 
the zero solution, so that (fj is non zero for even j only. Note that the 
inhomogeneous part hj of the j th order magnetic field amplitude Hj does not 
vanish for odd j. 

We study the long-time propagation by considering the unbounded or 
secular solutions, and a multi-time expansion. Therefore we introduce a 
sequence of slower and slower temporal variables T\ = r, r 2 , r 3 ,... defined 
by Tj = s 2 i +1 t. The propagation is governed by all equations of the KdV 
Hierarchy. In particular, the equation giving the evolution of the leading 
term <p 2 with regard to the first higher order time variable r 2 is derived as 
follows (more detail is given in jT2]). <£> 4 is the amplitude of the first correction 
to the main term whose amplitude is (f2- The equation that determines its 
evolution can be written in the form 

d T1 p A + qdt= (^2^4) + rd'lip A = -<9 T2 y? 2 - r 2 9f p a + 2 , (14) 
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where 2 refers to an expression depending on <p 2 without linear term, and r 2 
is a real coefficient. Some functions <^ 4 solutions of (fT4*|) are secular, i.e. grow 
linearly with the time T\. Consider values of the time variable t about 1/e 5 . 
Then the the time variable t\ = e 3 t take values about 1/e 2 , and the secular 
term in <^ 4 becomes of order e 2 instead of e 4 , due to the factor n oc 1/e 2 . For 
times with this order of magnitude, this correction term must be taken into 
account in the expression of the main amplitude (p 2 . In order to incorporate 
the correction into the evolution of the main amplitude ip 2 with regard to the 
second order time variable t 2 , we impose some condition on the rhs member 
of equation (JHJ), so that <y9 4 remains bounded (or more exactly sublinear). 
The condition to be satisfied is thus that the equation (J 14)) does not admit 
any secular solution. Through an explicit computation in the case where ip 2 is 
the one-soliton solution of KdV, Kodama and Taniuti JE] have noticed that 
the secular-producing terms are the terms linear with regard to the solution 
of lowest order (p 2 . The secular solutions ^4 will vanish thus if the linear 
terms vanish from the rhs member of the equation (|14|). To achieve this, we 
impose that (p 2 satisfies some partial differential equation such that 

d T2 <p 2 = -r 2 dl<p 2 + 2 . (15) 

We still need to determine the nonlinear terms of equation (|15jh represented 
by 2 . They are not free but imposed by the compatibility condition between 
the KdV equation © and the equation (|T5jl. which is the Schwartz condition: 
d Tl d T2 ip 2 = d T2 d T1 (p 2 . Kraenkel, Manna, and Pereira [IS] have conjectured and 
checked on many examples that the only equation that possesses the same 
homogeneity properties as the rhs member of (|14|). and that satisfies this 
condition, is the second equation of what is called the KdV Hierarchy. 
The KdV Hierarchy is the following family of equations ^7] : 

&T n v = dx£ n v {n integer), (16) 

where £ is a recurrence operator, defined by 

C = -\d 2 x -v + \ f X dX{d x v). (17) 

For n — 1, it is the KdV equation, with a normalization that differs from 
that of P (<? = §, r = A We identify both using the relations 

v = -fcp 2 , X = £ and T 1 = 4m. (18) 
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For n = 2, the equation of the Hierarchy writes as 

9t 2 v = \-d x v + -(d x v)d 2 x v + \vd\v + ^-v 2 d x v. (19) 
10 4 o o 

An important property is the existence of the r Hirota function [TZj, that is 
a function of all variables (X, Ti, T 2 , . . .), related to i> by 

T\, T 2 , . . .) = 2<9| In r(X, T\, T 2 , . . .) 

(avoid any confusion between the r Hirota function and the time variables 
Tj). The existence of r ensures that a solution v of the system yielded by 
all equations of the Hierarchy exists, and thus that the Schwartz condition 
is satisfied at any order. After an adequate choice of the proportionality 
constant that connects the time variables of order 2, the variable t 2 of our 
expansion and the variable T 2 of the Hierarchy, that we write as 

T 2 = -16r 2 r 2 , (20) 

the evolution equation to be satisfied by (p 2 is 

j^r d T2V2 = <%£V2. (21) 

This way, the linear terms have been removed from the equation (|15j). It 
remains to justify that this procedure, that removes all linear terms from the 
rhs member of the linearized KdV equation, assuming it polynomial with 
regard to the solution of KdV, ensures that the solution of the linearized 
equation is bounded [TH]. The KdV equation admits an infinite sequence of 
conserved densities we denote by Aj, an expression of which can be found 
in PHI- It has been proven in ^H] that the secular-producing terms are 
the terms proportional to d^Aj. Further the relations existing between the 
conserved densities Aj, and the recurrence operator C, defined by (|I7jL that 
allows to write the Hierarchy, allow to show that the procedure by Kraenkel 
et al, initially intended to remove the linear terms, exactly removes all these 
secular-producing terms. 

Otherwise, the rhs member of the linearized KdV equation that governs 
the evolution of <^6 involves (p&, solution of ((Hj). It is thus necessary to see 
wether, when a solution of the linearized KdV equation itself is used in the rhs 
member, which part of it is secular-producing, and which part is not. This is 



7 



not too difficult. Indeed, this solution is given by its expansion on the basis 
of the squares $fc of the Jost functions related to the KdV equation [TBJ EE] , 
and we have characterized the fact that a source term is secular-producing 
or not by some criterion, that involves the coefficients of this expansion and 
their t-dependency. It remains a last point to be studied: the dependency 
of the higher order terms with regard to the higher order times. We shown 
that it is governed by a linearized KdV Hierarchy [T3] . Finally, we have been 
able to justify that the higher order terms are not secular-producing, and to 
prove that the formal expansion contains bounded terms only. 



3 Time scales 

The generalization of the above procedure to an arbitrary order n ^ 2 yields 
the equation 

(Z^T^ = (22) 

which governs the evolution of the main amplitude ^2 with regard to the 
higher order time variable r n . £ is defined by the above formula (|17jl. The 
scaling coefficient r n is defined by 7*1 = r and the recurrence formula 



n-l 

■ n+1 = m^)i^n-i, k)l[(-r^. (23) 

7=1 

( a i)i<j«»-i. 
(E"=i 2ja j )+k=2n+3 

The sequence of time variables T\, r 2 , t 3) ... involved by the multiple 
time formalism are thus affected by the sequence of scaling coefficients r\, 
T2, rs,... . The equations of the KdV Hierarchy are 'universal', not specific 
to the physical situation considered. The time scaling coefficients contain 
thus most physical data about the time evolution of the wave. Further, they 
are of interest regarding the convergence of the asymptotic series. They are 
computed using recurrence formula ((2HJ), together with the results of [12 j 
listed in the appendix. The first coefficients read as follows: 



r 2 



8 + (4a - 13)7 + (3a + 6)7 2 



8(1 + a) 2 m^/i 



(24) 

+a{Aa - 10)7 3 + 4a(l - a)^ A 
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4T/13 



r 3 = ian ?\s 62 40 + ( 32a - 88 h + ( 8 « 2 + 12 « + 67 h 2 
16(1 + a.) 6 rrqix z L 

+ (68a 2 - 132a - 20)7 3 + (40a 3 - 143a 2 + 122a + 2) 7 4 
-a(32a 2 - 70a + 32) 7 5 + a 2 (16a 2 - 56a + 6)7 6 

-8a 2 (4a 2 - 8a + 1)7 7 + 16a 3 (a - 1)7 8 



(25) 



in which 
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li = 1 + «7 



m t = msin^. 



(26) 



The expressions of the higher order coefficients can be obtained in the same 
way, but are too complicated to be written down here; numerical computation 
is more convenient. 

The coefficients r\ = r, r 2 ,... up to r 5 are plotted on figured against 




Figure 1: Plot of the five first time scaling coefficients r±, r$ against the angle 9 
between the propagation direction and the exterior field. The rescaled magnetic induction 
is m = 1, and the parameter determining the strength of the exterior field is a = 0.5 . 
Dotted line: r\, solid line: r%<, large dashing: 7-3, dashed-dotted line: r±, short dashing: r§. 
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the value of the angle 9 between the propagation direction and the external 
field, for a given value of the parameter a that determines the magnitude of 
this field. Notice the annulation and sign change of the coefficients and 
rs, about 0.41 and 0.48 radians respectively. This marks a change in the 
behaviour of the corresponding corrections. When is zero, <£>2 is constant 
with regard to r 4 , thus the 3 rd order correction is in fact valid at order 4, 
regarding its time dependency. 

It is seen that the r n take very small values when 9 is close to tt/2, and 
very large values when 9 is small. In the limiting case where the propagation 
direction is orthogonal to the external field (9 = tt/2), the velocity V is 1, 
thus 7 = 0, and the coefficients q and r vanish, so that the KdV equation 
© is replaced by 

d T ip = 0. (27) 

Thus if is constant with time at fist order, which means a priori that the 
wave evolves much slower than in the general case, at least for an order 
of magnitude. Recall that this order of magnitude is determined by the 
perturbative parameter e, related to the wave amplitude and typical length. 
The wave propagates without deformation, to within a quantity of higher 
order in e, up to times about T/e 3 , instead of times about T/e, as usual in 
the long wave approximation. The higher order equations simplified by this 
trivial time evolution of the main term, yield also an approximation valid 
up to T/e 3 , for some finite T. Let us precise the influence of the scaling 
coefficients on the time validity range of the higher order approximations. 
We denote by Lq some typical length of the wave. The dimensionless space 
variable is £,/Lq = e(x — Vt)/Lo- The reference length for x is chosen with 
the order of magnitude of eL Q , in such a way that, as x takes values as large 
as 1/e with respect to this reference length, £ is about L . V is close to 
1, thus taking the same value as a reference time (recall that t has already 
been rescaled into ct) is coherent with the asymptotic expansion. The higher 
order time variables adapted to the expansion are, rather than the r n , the 
variables T n of the KdV Hierarchy (|2*2*|) . written under its normalized form 

dr n v = d^C n v (n integer). (28) 

The differential recurrence operator C is as in (fTTj) . with v = • The 

variable T n reads then: 

T n = -{-A) n r n r n = -(-A) n r n e 2n+ h (29) 
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T n must be about 1 for large values of t. This necessitates a smaller value 
of e when the coefficient r n , or rather (— 4) n r n , is large. For each n, e must 
be compared to the reference value e n defined as follows: e n is the value of e 
in (j2Hj) such that, when T = et is equal to L Q , \T n \ has the same value. It 
yields 



The approximation involving the first n time variables (for all written terms) 
is valid for t T/e 2n+1 , for some finite T with the order of magnitude of 
unity in the initial unit (eLqv). Taking the scaling coefficients into account, 
the approximation will rather be valid for \T n \ ^ T, that is: 



Notice that it is in fact necessary that \T p \ ^ T for all p ^ n, what implies 
some conditions on the variations of e p with relation to p. According to (|31jl , 
when e n takes large values, the propagation can be described over a long 
distance even if the order n is relatively low and the value of the perturbative 
parameter e close to 1. The higher order time variables make sense only if 
e is smaller than the e n , and long time propagation can be described only if 
the ratio e/e n is very small. These conditions will be hardly fulfilled when 
e n becomes small. 

The five first Eyi £1X6 plotted on figure El against the angle 9, and on figure 
H3 against the ratio a that determines the magnitude of the external field. If 
the extrapolation of the few computed terms is valid, the sequence e n seems 
to be bounded with regard to n, although its terms grow up as 9 tends to 7r/2. 
Further, this bound is not excessively small when 9 is not smaller than 10° or 
15°. For smaller values of 9, the e n are so small that the KdV approximation 
can be valid only for excessively low intensities, and the higher orders will 
never appear. 

When 9 approaches n/2, the e n become large. Then the approximation 
yielded by the KdV equations is valid during a very long time. The pulse 
behaviour will be correctly described by them even if the small perturbative 
parameter e takes values rather close to 1. At the limit 9 = 7r/2, the modu- 
lation described by the KdV equation itself arises only at a very slow rate. 
A typical dependency of the e n with regard to the strength of the external 
field is shown on figure H3 The e n grow slowly with a. Thus a strong ex- 
ternal field enhances the validity of the KdV approximation, and increases 




(31) 
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Figure 2: Logarithmic plot of the five first e n , reference values for the perturbative 
parameter s, built from the time scaling coefficients r n , against the angle 9 between the 
propagation direction and the exterior field. The value of the constants and the legend 
are the same as in figure ^ 



duration along which it can be expected to describe the physics. However 
this effect is much weaker than the dependency with regard to the direction 
of the external field and the angle 9. 

4 The soliton of the hierarchy 

The time scaling coefficients studied in the previous paragraph have given 
an insight into the convergence of the perturbative expansion as an asymp- 
totic behaviour for small values of the perturbative parameter e, for a fixed 
number n of corrective terms. We are also able to get some insight into the 
convergence of the series when n tends to infinity and e is fixed, through the 
computation of the one soliton solution of the complete KdV hierarchy. As 
mentioned above, all equation of the KdV hierarchy are compatible together, 
in the sense that for a given initial data, a function v(X, T\, T2, T3, ■ ■ •) sat- 
isfying equation (|TT)j) for any value of n can be found. This solution can 
be found using the inverse scattering transform (1ST) method, at least in 
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O.lL 

Figure 3: Logarithmic plot of ei, £5, against the parameter a that determines the 
strength of the exterior field. The rescaled magnetic induction is to = 1, and the angle 
between the propagation direction and the exterior field is 9 = n/A . The legend is the 
same as in figure ^ 



principle. Indeed, all equations of the hierarchy are completely integrable by 
means of the 1ST method. Furthermore, they can all be described in the 1ST 
formalism using the same spectral problem (^H], P- 96), which ensures their 
compatibility. The scattering data (R + (k),D + j,kj) (see [TH], p. 141 sq., for 
the precise definition of these quantities) are defined in the same way for all 
equations, only their time evolution differ for each time variable T n . These 
time evolution is given by ([19J, p. 149) 

R + (k,T n ) = i? + (M)e^ (fc)T ", (32) 
D +J (T n ) = D +ij (0)e n ^ Tn , (33) 
k,{T n ) = kj(0). (34) 

The index n refers to the n th equation of the hierarchy. The evolution factors 
are fl n j = fl n (kj), and fl n {k) = —iujj(2k), where uij{k) is the dispersion 
relation of the n th equation of the hierarchy linearized. Its seen from relation 
(|3~3jl that the discrete spectrum (kj) is constant with regard to any of the 
time variables T„ . Therefore the number of solitons and their characteristics 
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are not modified by the higher order time evolution. The evolution of the 
spectral data with regard to all the higher order time variables can then be 
written as a single exponential factor for each spectral component, 



R+{k,T 1: T 2 ,---) =i2 + (fc,0,0,---)exp ^fi n (fc)T n . (35) 



v n=l 



From the expression ()16II17|) of the equations of the hierarchy, we find that 

— k 2n+1 

^n(k) = (36) 

For a value of the spectral parameter k belonging to the discrete spectrum, 
k = kj = iKj with Kj real, we get 

n nJ = 2{-l) n+1 K f +1 . (37) 

Using the definition ()29j) of the time variable Tj , we get the following expres- 
sion of the complete time evolution factor: 

oo oo 

tt nJ T n = Qjt, with Qj = Y (2en j ) 2n+l r n . (38) 

?1=1 71=1 



Obviously formula ()38|) is valid only if the power series converges. Notice that 
the coefficients of the latter are the time scaling coefficients r n . For a one- 
soliton solution, the above formulas show that the introduction of a sequence 
of higher order time variables and of all equations of the KdV hierarchy 
yield nothing but a renormalization of the soliton speed. This result can be 
also found by direct computation as follows. By definition, the one-soliton 
solution propagates without deformation, at least with regard to the first time 
variable 7\. It can thus be written under the form v = v(X + ATi). Then 
using the KdV equation, i.e. equation (|T6|) with n — 1, we see that v is an 
eigenvector of the recurrence operator L defined by fll7j) . with the eigenvalue 
A. We deduce easily the T n -dependency of v, it is given by v — v(X + X n T n ). 
We find this way the expression of the one-soliton solution of the complete 
hierarchy: 

v = 2b 2 sech 2 b (x + jr{-b 2 ) n T n ) , (39) 
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using the normalized variables. In the case of magnetic solitons, it can be 
written using the physical variables as 



1 0r -> 

H w = —hi(3 2 sech 2 (3{x - Vt) , (40) 



where 

V = V 4- 



V = V + ^A n p 2n r n . (41) 

n=l 



V is the velocity given by ©, hi the polarization vector given by fTOj) . The 
wave magnetic field if™ is related to the previously defined field components 
through 

H = H + e 2 H 2 + ---~H + H W . (42) 

The dimensional soliton parameter (3 is related to the normalized soliton 
parameter b through (3 = eb. Computation of the one-soliton from the 1ST 
formalism allows to identify the soliton parameter b to the single discrete 
eigenvalue K\. This way we check that the relative soliton velocity (V — V) 
given by (j4l} is equal to f2i/(2/3), using the expression (J38|) of the evolution 
factor Q,i. 

The soliton speed is thus given by a power series of the soliton parame- 
ter (3, whose coefficients are essentially the time scaling coefficients (r n ) n ^i. 
Obviously if this series diverges so does the whole perturbative scheme. Re- 
ciprocally, the convergence of the series defining the velocity should favour 
that of the perturbative scheme, although the latter is by no means proven. 
Writing the power series which defines V as 

n=l Vfc ™/ 

we see that it converges when 

P < I3 M = lim inf e n (44) 

n >oo 

and diverges for larger values of the soliton parameter f3. Therefore the 
limit of the sequence e n for large n gives us a maximal value (3m of the 
soliton parameter (3, above which we know that the perturbative scheme 
does not converge when increasing the number of terms. Physically, this lack 
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of convergence means that the KdV soliton will be destroyed by some effects 
which cannot be taken into account using the perturbative approach. 

For values of the soliton parameter (3 below the limit (3 m, we get a renor- 
malized soliton speed, a priori valid for an infinite propagation time. The 
boundness of all terms in the perturbative scheme proves that for a given 
propagation time and a given number of terms, this soliton gives a good 
approximation of the real impulsion for small enough values of e, i.e. of (3. 
We can reasonably conjecture that small enough can be understood here as 
less than the limiting value (3m of (3. Physically it means that magnetic KdV 
solitons with parameter smaller than (3m should conserve their properties 
during a long propagation time. 

According to figures 121 and El the e n , thus also their limit (3m, depend on 
the physical parameters, and specially on the angle 9 between the propagation 
direction and the applied field. An example of computation showing the 
convergence of the velocity series is drawn on figure 0] as a function of this 
angle. It is seen that, for 9 close to tt/2, the first approximation (KdV) gives 
almost the exact speed, while for small angles the series diverges. We denote 
by 9 M the value of 9 for which (3 M is equal to the fixed value of (3. When 
9 < 9m, the series does not converge, and the whole perturbative approach is 
not valid. To compute 9m for the figure El we have approximated (3m by £5. 
When 9 > 9m, if we consider only the soliton speed, the KdV approximation 
will correctly describe the wave evolution. More precisely, the KdV equation 
itself will give an acceptable description above some value 9 t of the angle 9, 
while this first order approximation needs to be corrected by higher order 
terms below 9 t (notice that the threshold value 9m is precisely defined, while 
9t is only an order of magnitude depending on the accuracy required). It is 
reasonable to think that the same kind of conclusion holds in a more general 
situation, involving several solitons and radiation. 

5 Conclusion 

The multiple time formalism has been applied to the study of the propaga- 
tion of KdV solitons in ferromagnetic media. According to this formalism, 
the dependency of the higher order terms with respect to the first order 
time variable is given by linearized KdV equations, while the dependency 
of the main term with regard to the higher order time variables is governed 
by all equations of the KdV Hierarchy. The latter are determined by the 
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Figure 4: Plot of the five first approximate values V n — V + Y^>=i (P/ £ n) 2n of the soliton 
velocity V, against the angle 8 between the propagation direction and the exterior field. 
The rescaled magnetic induction is m = 1, and the parameter determining the strength of 
the exterior field is a = 0.5 . Solid line: n — 0, dotted line: n = 1, large dashing: n = 2, 
dashed-dotted line: n = 3, short dashing: n — 4. For a soliton parameter (3 — 0.1 (a), 
p = 1.5(b). 



requirement that the linear terms in the rhs of the linearized KdV equation 
vanish. This yields scaling coefficients for the higher order time variables of 
the KdV Hierarchy, which contain most physical information concerning the 
wave evolution. Explicit computation of these scaling coefficients shows in 
particular that the approximation yielded by the KdV model gives a good 
account for the physical behavior of the wave during long propagation times 
when the angle between the propagation direction and the external field is 
large enough. The time during which the pulse is correctly described by 
the KdV equation falls to zero when they are parallel. Mathematically, the 
perturbative parameter e is infinitely small, while it takes a finite value in 
a physical situation. The approximation is valid only if this finite value is 
small enough. The corresponding range of the perturbative parameter e is 
usually determined in a rather empirical way. The present study gives some 
theoretical insight into this question, through a physical interpretation of the 
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time scaling coefficients. 

The one-soliton solution of the complete KdV hierarchy has been written 
down as a function of the physical parameters. The soliton velocity writes 
as a power series of the soliton parameter, involving the sequence of the 
time scaling coefficients. We get a maximum value of the soliton parameter, 
above which the perturbative series diverges. Then the KdV approximation, 
even with corrective terms, does not describe correctly the physics. If the 
soliton parameter is below the threshold, the long-distance effect of the higher 
order corrections is only a modification of the soliton speed, and the physical 
system behaves qualitatively as the KdV model. It has been observed that the 
validity domain of KdV-type asymptotics is often much larger than predicted 
by the mathematical analysis. The above conclusions can partially explain 
this observation: the KdV-type behaviour is qualitatively correct in the whole 
validity domain of the infinite KdV hierarchy expansion, which is expected 
to be much larger. 



18 



Appendix 



We list in this appendix the formulas needed for the computation of the 
time scaling coefficients r n . These formulas are proven in [T2j. r n is given by 
equation with 



E{(aj) m , k ) = -^ [Vm ■ rh ((o^i, k-1) 
- J^m • m ((atj - .,. k-1) 



where 



m 



m x 
m t 




m {{ctj), k) is deduced according to 

/ e((aj) m ,k,l) 
u ((aj) m , k, Z) = h ((aj) m , k, I) 
\ m ((aj) m ,kj) 

from the following recurrence formulas. For all / ^ 1, 

u ((0), 0, Z) = 

For all k and 1^1, 

u ((Q),k,l) = S(Vm ((0),fc- 1,Z)). 
For all 7^ (0) and 1^1, 

u ((a^i, 0, Z) = 

((a,- - 0, Z)) 



For all (c^)j>i ^ (0), fc, Z ^ 0, 

S (Vm ({aij)fei,k- 1,0) 

+ ((oj - <y w )j>i, 0, z)) 



-5" (m ((a,- - fc - 1, Z)) 



(45) 



(46) 



(47) 

(48) 
(49) 

(50) 



(51) 
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In 



, Si j is the Kronecker symbol, V is given by (0). $ is denned by 



$(it) = S(arh A rhj(u)) + Uj(u) 
where S is the 9x3 matrix 

-hR 



S = TL- 1 with T = 

I is the three-dimensional unity matrix, and 

1 



I 

-r 



L- 1 



\im x m t 

Ui is the linear operator in M 9 defined by 


















m t 


m x 









with 



E 

»/(| # 

M 



E 

)Tn( I # 



.1/ 







-Ra; is the 3x3 matrix 




0-1 
1 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



The first term of the sequence u ((o£j)j^i, k, I) is given by u\ = Thi, where 
hi is the polarization vector defined by (fTUj) . which also reads 



fxm x 
1 + a)m t 




(58) 
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The quantity A in (J45|) is given by 

A = Vrh ■ $ m (tti), (59) 
where $ m is the m-component of $ defined by ([52)1. according to 

(60) 




We use the shortcuts 




m x = mcosB, (61) 



and 7, fi, m t given by (}2l)j) . The expression (0) of the velocity yields the 
relation 

fiml + 7(1 + a)m t 2 = 0, (62) 
which is useful to simplify the expressions. 
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